Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 2 February 2008 (MN KTeX style file v2.2) 



Constraints on dark energy models from galaxy clusters with 
multiple arcs 

Massimo Meneghetti^, Bhuvnesh Jain^, Matthias Bartelmann^ and Klaus Dolag-^ 

^ITA, Universitdt Heidelberg, Tiergartenstr. 15, D-69121 Heidelberg 

^Department of Physics, University of Pennsylvania, Philadelphia, PA 19101. USA 

^ Dipartimento di Astronomia, Universita di Padova, Vicolo deU'Osservatorio 2, 1-35120 Pudova 



2 February 2008 



ABSTRACT 

We make an exploratory study of how well dark energy models can be constrained using 
lensed arcs at different redshifts behind cluster lenses. Arcs trace the critical curves of clus- 
ters, and the growth of critical curves with source redshift is sensitive to the dark-energy 
equation of state. Using analytical models and numerically simulated clusters, we explore the 
key factors involved in using cluster arcs as a probe of dark energy. We quantify the sensitiv- 
ity to lens mass, concentration and ellipticity with analytical models that include the effects 
of dark energy on halo structure. We show with simple examples how degeneracies between 
mass models and cosmography may be broken using arcs at multiple redshifts or additional 
constraints on the lens density profile. However we conclude that the requirements on the 
data are so stringent that it is very unlikely that robust constraints can be obtained from in- 
dividual clusters. We argue that surveys of clusters, analyzed in conjunction with numerical 
simulations, are a more promising prospect for arc-cosmography. 

We use such numerically simulated clusters to estimate how large a sample of clus- 
ters/arcs could provide interesting constraints on dark energy models. We focus on the scatter 
produced by differences in the mass distribution of individual clusters. We find from our sam- 
ple of simulated clusters that at least 1000 pairs of arcs are needed to obtain constraints if the 
mass distribution of individual clusters is taken to be undetermined. We discuss several un- 
solved problems that need study to fully develop this method for precision studies with future 
surveys. 
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1 INTRODUCTION 

The compelling observational evidence for an accelerating expan- 
sion of the universe has led to much recent work on models for 
dark energy that would drive this expansion (see Peebles & Ratra 
[2002; Padmanabhan 2003 and Carroll 2003 for reviews). Unless the 
dark energy is a cosmological constant, it is likely that it evolves 
with time. The best observational probes of dark energy are there- 
fore ones that measure the geometry of the universe at different 
redshifts, allowing for a mapping of the evolution of dark en- 
ergy. So far Type-la supernovae, which measure the luminosity 
distance over a range of redshifts, and the multipole orders of the 
acoustic peaks in the power spectrum of the cosmic microwave 
background , have been the only geometric probes of dar k en ergy 
jPerlmutter et al..l999.: .Riess et al... 1998; -Spergel et aL.26o3l : see 
iHuterer & Turneil200ll for a discussion of different probes). 

The multiple imaging of background galaxies into arcs by 
foreground galaxy clusters has been used to probe the mass dis- 
tribution of clusters. The position of the arcs on the sky depends 
on the mass enclosed within the lensing region of the cluster and 



also on ratios of angular diameter distances, which we will refer 
to as the geometric information probed by lensing. Observations of 
multiple arcs at different redshifts have been suggested as probes o f 
dark energy jU ink & Piercell998l:lGolse et all20 02VSereno'2002h. 
The idea is that the relative positions of arcs at different redshifts 
depend only weakly on the mass dist ribution and therefore probes 
the lensing geometry. Recently S oucail et alJi200 4) analysed Abell 
2218 using arcs at four different redshifts to constrain dark energy 
models, obtaining w< —0.85 at 68.3% confidence for constant w. 

A critical issue in the reliability of geometric measurements 
from cluster lensing is the sensitivity of the results to the mass 
distribution of the cluster. If the cluster mass were smoothly dis- 
tributed, then it is easy to see that multiple arcs are a good probe 
of geometry, but realistic clusters are likely to have substructure 
and ellipticity. Ideally one would like to be able to independently 
measure the mass distribution responsible for lensing, so that no 
theoretical assumptions need be made in inferring geometric infor- 
mation (e.g. C hiba & Taka hashi 2002)- 

In this paper we use numerical simulations of lensing clusters 
to explore two approaches to constraining dark energy models with 
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cluster arcs. The first may be called the "golden lens" approach: 
we ask how external information about the cluster mass distribu- 
tion can allow dark energy constraints to be obtained from just one 
or a few lens systems. The second is a statistical approach: we ask 
how large a sample of clusters is needed to make an "ensemble av- 
eraged" measurement of critical curve versus redshift, which can 
then be compared to simulations. The two approaches seek differ- 
ent kinds of datasets, and the latter assumes that simulations (cali- 
brated to some extent from the data) represent lensing clusters ade- 
quately well at least for the giant arcs. 

We describe our parameterisation of dark energy models in 
Sect. 2. The lensing formalism is presented in Sect. 3, and the 
analytic dependence of critical curves on cosmology is shown. In 
Sect. 4, we extend the study to numerical cluster models, and we 
suimnarise in Sect. 5. 



Table 1. Parameters characterising the cosmological models 



2 DARK-ENERGY MODELS 

We work with the metric 



where we have used the comoving coordinate % and fl(x) 



(1) 
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z)^ is the scale factor as a function of conformal time x. We adopt 
units such that c = 1. The comoving angular diameter distance 
depends on the curvature: we assume a spatially flat universe so that 
r(x) = X- The density parameter SI has contributions from mass 
density flm or dark energy density iJje, so that SI = flm + flje. 

Dark energy models are often described in terms of the equa- 
tion of state p = wp, with = — 1 corresponding to a cosmological 
constant. The time dependence of w is commonly parameterised as 
w = H'o + vi'a(l — a) ( Linder.2003.) . For comparison with other work 
we will compare Wa to w' defined by w = wq + w'z- The Hubble pa- 
rameter H{a) is given by 

1 /2 

i/H=/fo[fim«"^+^idee"^-^° , (2) 

where //q is the Hubble parameter today. The comoving distance 
X(«) is 



X(«) 
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The above equation shows how the lensing observables depend on 
integrals over the dark-energy dependent expansion rate. 

We consider four different cosmological models in this study 
(see Table 1 and iDolag et aL.2004, for details). These can be de- 
scribed in terms of an effective values of wq and Wa- One is 
the standard cosmological-constant model with flat geometry and 
Slj^ = 0.7. The second has a constant ratio wo = —0.6 between 
pressure and energy density of the dark energy (DECDM). The re- 
maining two use different descriptions for the self-interaction po- 
tential of the dark-energy scal ar field. One has a pow er-law po- 
tential (Ratra-Peebles, RP; see IPeebles & R'atr3l2002h . the other 
has a powe r-law potential multiplied by an exponential (SUGRA; 
iBrax & M artin 2000). The RP and SUGRA models have the same 
values of w at the present time, wq — 0.8, but very different evolu- 
tion histories. See lDolag et alJ^2004^ for other details of the models 
and their implementation in the simulations used here. 
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3 ANALYTIC MODELS 
3.1 General properties 

We adopt the density profile proposed by iNavarro et alj ^993) 
(hereafter NFW) for modelling cluster lenses, which was found 
to fit numerically simulated galaxy clusters well. We investigate 
here how strong lensing by analytical NFW halos changes when 
the dark-energy equation of state is changed, and postpone effects 
of the substructure and asymmetry of realistic halos to Sect. 5. The 
profile is given by 



rA,)(l + rA,)2 



(4) 



where Ps and rs are characteristic density and distance scales, 
respectively. These two parameters are strongly correlated 
jNavarro et al.ll997^ . 

The concentration is defined as c = r2Qo/r^, where r2()o en- 
closes a mean halo density of 200 times the critical cosmic density. 
Numerical simulations show that c depends on the virial mass M of 
the halo, which can thus be used as the only free parameter. Sev- 
eral algorithms have been developed for relating c to M. In this 
work, we adopt that pr oposed bv, E ke et al., 12001) because, as re- 
cently demonstrated bv lDolag et alj i20o'4l) . it performs very well 
also in dark-energy cosmologies. The concentration also depends 
on the cosmological model, implying that the lensing properties 
of haloes with identical mass are different in different cosmologi- 
cal models if they are modelled with the NFW pro file. This is an 
important difference of our work to e arlier studies iLink & Piercd 
ll998tlGolse et all2002tls'erenj2002h . 

The lensing properties of the NFW p rofile have been 
widely investigated in the past (see e.g . Bartelmann 199g; 
[ Wrigh t & Brain er j20()(ilLi & Ostrikedl2002:iPerrotta et al.„20ol 
[Meneghetti et alT l2003fr thus we summarise them only briefly. 
Lensing is fully described by the lensing potential \|/. For axially 
symmetric models, computing the lensing potential reduces to a 
one-dimensional problem. We define the optical axis as the straight 
line passing through the observer and the lens centre and intro- 
duce the physical distances perpendicular to the optical axis on 
the lens and source planes, ^ and T], respectively. We then choose 
Ts as a length scale in the lens plane and define the dimension- 
less distance x = ^/r, from the lens centre. By projecting to the 
source distance, we define a corresponding dimensionless distance 
y = (r|/rs)(£)i/Ds) from the optical axis in the source plane. 

Defining = p,r,I.-,\ where = [cV(47tG)] [£)s/(OiAs)] 
is the critical surface mass density for strong lensing and Djs is the 
angular diameter distance between the lens and the source planes, 
the lensing potential can be written as 



where 



(5) 
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Axially symmetric models are g enerally inappropriate fo r de- 
scribing lensing by galaxy clusters fMeneshetti et al. 2003) be- 
cause their typically high degree of asymmetry and substructure 
changes their lensing properties qualitatively and substantially. The 
tidal (shear) field produced by an asymmetric mass distribution can 
be partially mimicked by incl uding ellipticity in to the model. We 
adopt the model proposed by iMeneghetti et alj j2003), who ob- 
tained an elliptical generalisation of the NFW axially symmetric 
lensing potential by substituting in Eq. Js} 
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+ x2(l-e) 
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where xi and X2 are the two Cartesian components of x, x^ = x^ + 
x^. We define the ellipticity as e = l—b/a, where a and b are the 
major and minor axes of the ellipse, respectively. 

The lensing potential implies the Jacobian matrix of the lens 
mapping, 



d'^Mf{x) 
dxjdxj 



(8) 



whose determinant vanishes on the critical curves of the lens. In 
particular, the tangential critical curve is located where the tangen- 
tial eigenvalue of the Jacobian matrix vanishes, 

:^ = 1-K-Y=0, (9) 

where K and Y= (Yi i Y2) are the lens convergence and shear, respec- 
tively. They can be written in terms of the lensing potential as 

1 



K(X) 

and 
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llix) = ¥12 =¥21, 
where we abbreviate 
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3.2 Sizes of critical curves and cosmology 

Tangential arcs form near the tangential critical curves of a lensing 
cluster. We assume in the following that there is a correspondence 
between the tangential critical curves and the observed position of 
tangential arcs. In practice, tangential arcs with moderate length- 
to-width ratio may form quite far from tangential critical curves if 
the cluster is embedded into a strong shear field. To reduce this un- 
certainty, arcs with a large tangential-to-radial magnification ratio 
must be used. 

The problem is illustrated in Fig. Q where we show the po- 
sition of all images with length-to-width ratio > 5 (left panel) and 
> 10 (right panel), found in a ray-tracing simulation using a numer- 
ical cluster model as a lens. A large number of elliptical sources 
were distributed in the region containing the lens caustics in or- 
der to obtain a large number of arcs. In this particular case, the 
tangential critical curve appears very elongated, showing that the 
cluster mass distribution is far from axially symmetric. From the 
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Figure 1. Position of gravitational arcs with length-to-width ratio larger 
than 5 (left panel) and 10 (right panel) in a lensing simulation; also shown 
are the critical curves of the lensing cluster 



first panel it is clear that arcs with substantial tangential distortion 
form in a wide region surrounding the critical curves. For example, 
using arcs with a length-to-width ratio > 5, the spread in arc posi- 
tions is ~ 10"; selecting arcs with high length-to- width ratio (> 10) 
reduces the spread to a few arc seconds. 

It is important to note that arcs with large length-to-width ra- 
tios form along those portions of the critical curves which are at 
the largest distances from the cluster centre. In comparing theory 
to observations, one possibility is to compare in detail the arcs in a 
simulated sample with observed arcs. However to avoid the compu- 
tational expense of lensing source galaxies through the numerical 
clusters discussed in the next section, we will use an average size 
of the critical curve weighted with the inverse radial magnification 
as a proxy for the location of large arcs. Hence we do the same 
here by assigning to the i-th point on the tangential critical curve a 
weight: 



Pi 



|2[l-K(f,0]| 



(14) 



The size of the critical curve is then measured as the weighted av- 
erage distance of the critical points from the lens centre. 



Li Pi 



(15) 



The weight factor that we have used typically underestimates the 
location of giant arcs. However we use it for simplicity since we 
do not expect it to significantly change the relative sizes of critical 
curves at different redshifts. DalaLet_al. (2003) give a more detailed 
discussion of the location and orientation of giant arcs in terms of 
the profiles of K and y. 

We assume here that the cluster lens is at redshift z\ ~ 0.6, and 
shift the source plane from Zs = 1 to Zs = 5. The virial mass of the 
lens is 7 X lO'"* h^^ Mq. Using Eq. <15> . we measure the size of the 
critical curves produced by our analytic model for different source 
redshifts. This is repeated for all cosmological models we consider. 
For each of them, we finally estimate the growth rate of the lens 
critical curves. 

We show in Fig.|2|how the critical curve grows with source 
redshift. The curves are normalised to the size of the critical curve 
for sources at Zs = I. Thin and thick lines show the results for el- 
lipticities e = and e = 0.3, respectively. The growth of the critical 
curve with redshift is larger in the ACDM model than other cos- 
mologies with dark energy. In this case, for the axially symmetric 
lens, the growth between Zs = 1 and Zs = 2 is roughly by a factor of 
4.5. Beyond Zs = 2, the growth slows down; for sources at Zs = 5, 
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Figure 2. Sizes of the critical curves for sources at Zs. normalised to their 
value at z., = 1. Results are shown for four cosmological models. The lens 
mass is 7 x lO'* ft^'M,., and the lens redshift is z\ = 0.6. The lens is mod- 
elled using both an axially symmetric (thin curves) and a pseudo-elliptical 
(heavy curves) model with NFW density profile. In the elliptical case, the 
ellipticity of the iso-potential contours is e = 0.3. 



critical curves are larger by a factor ~ 7 than for sources at Zs = 1. 
Note that this increase is determined not just by the larger distances 
to higher redshifts, but also by the steepening of the NFW profile 
at large radii which are probed by the higher source redshifts. The 
model which deviates the most from ACDM is the SUGRA model, 
for which the critical curves for Zs = 2 and Zs = 5 are larger by fac- 
tors of ~ 3.6 and ~ 5.2, respectively, than for Zs = 1- The RP and 
the DECDM models fall between them. 

Adding ellipticity to the model reduces the growth rate of the 
critical curves for all cosmologies. For e = 0.3, the critical curves 
are larger by a factor of ~ 3.4 for Zs = 2 compared to Zs = 1 in 
the ACDM model and by a factor of ~ 2.8 in the SUGRA model. 
The dependence of the growth rate on the lens ellipticity will be 
discussed in detail in the following section. 



3.3 Dependence on lens mass, redshift and ellipticity 

As noted above, the amount by which the critical curves grow as 
a function of source redshift depends on the lens ellipticity. In this 
section, we discuss the dependence of the critical-curve growth rate 
on lens ellipticity, mass and redshift. 

To explain the dependence on these lens properties, we must 
consider where the critical curves form. In the case of an axially 
symmetric model, they occur where the radial function F{r,z) ~ 
K(r,z) +'Y{r,z) is unity. For the NFW density profile, F{r,z) can 
locally be parametrised as 



where the |3(r, z) is the logarithmic slope at radius r, 
aiogF(r,z) 



3logr 



(16) 
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Figure 3. Dependence of the critical-curve position on lens redshift, mass 
and ellipticity. Critical curves arise where log(f ) = 0, as indicated by the 
horizontal dotted lines. Shown is log(f ) as a function of log(/') for two lens 
redshifts (top left panel), for two lens masses (top right panel) and for two 
lens ellipticities (bottom left panel), hi the bottom right panel we show the 
ratio of the critical-curve size for sources at redshift Zs = 2 and Zs = 1 as a 
function of the logarithmic slope p measured on the Zs = 1 critical curve. 



which increases with radius. 

The dependence of F on the lens redshift, mass and ellipticity 
is shown in the first three panels of Fig. |3| In each plot we show 
log(F) as a function of log(r). Where the ellipticity is varied, we 
plot F along the major axis of the ellipse. When not explicitely 
mentioned, the lens model has a mass ofM = 7x lO'4/,-lM0,a 
redshift of zi = 0.6, and an ellipticity of e = 0. 

Keeping the source redshift fixed at Zs = 1, changing the lens 
redshift, mass or ellipticity causes the curve to be shifted up and 
down, or left and right, implying that the critical curve forms at 
radii characterised by different values of (3. In particular, when 
moving the lens closer to the sources or the observer, or decreasing 
its mass and its ellipticity, F reaches unity at larger p. 

The relative shift of the critical-curve position as a function of 
the source redshift depends on the local slope of log(F) where the 
critical curve forms. Since the tangential critical curve for sources 
at redshift z occurs where F{r,z) = 1, the critical curve for sources 
at redshift z' = z + dz lies where 



f(^r-\^')=fl^r;z') + 



dFir,z'] 



dr= 1, 



from which we obtain 

SFir,z') 



Ar=l-F{r,z')=AF{r,z'). 



dr 
Since 

dF _ d\nF F 
dr dlnr r ^F 
the relative shift of the critical curve is given by 
Ar _ 1 AF 



(18) 



(19) 



(20) 



(21) 



When log(F) is flatter, a larger relative shift of the critical curve is 
obtained. This is shown in the last panel of Fig.|3| where the ratio 
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Figure 4. Growth of the critical curve size as a function of the lens redshift, mass and ellipticity. The ratio between the sizes of the critical curves at redshifts 
Zs = 2 and = 1 is shown. 



between the sizes of the critical curves for sources at redshift Zs = 2 
and Zs = 1 is plotted as a function of P measured at the position of 
the critical curve for sources at Zs = 1- Here we have used as a 
lens an axially symmetric model of mass M = 7 x lO'^A^' Mq at 
redshift zi = 0.6. 

Since, as shown earlier, the value of P on the critical curve 
depends on the lens redshift, mass and ellipticity, we expect a dif- 
ferent growth rate of the critical curves for lenses with different 
values of these parameters. 

This is shown in Fig. |4] We plot the ratio of the sizes of the 
critical curves for sources at Zs = 2 and Zs = 1 as a function of 
lens redshift, mass and ellipticity. These are shown for the four 
dark energy models as indicated. As before, where not explicitly 
indicated, the lens model has mass M = 7 X 10'"* h-^MQy, redshift 
zi = 0.6, and ellipticity e = 0. As expected, the growth rate is larger 
for lenses at higher redshift and with lower mass and ellipticity. A 
stronger dependence is found on redshift and ellipticity, while the 
dependence on mass is much weaker. 

We notice that the sensitivity of the displacement of the crit- 
ical curves to these lens properties changes as we vary the equa- 
tion of state of dark energy. This reflects the different evolution 
of halos of a given mass in d i fferent cosmologies . As shown by 
iBartelmann et a l.' l"2002 . '2003^: lDolag et al.l i2004h . the formation 
epoch of dark-matter halos with the same mass in these cosmolog- 
ical models is significantly different, leading to substantially dif- 
ferent concentration parameters of the halo density profiles. Halos 
tend to form earlier and thus have typically larger concentrations 
in the SUGRA and in the DECDM models compared to the RP 
and the ACDM models. Thus, for these models, the growth rate of 
the critical curves is less sensitive to changes in redshift, mass and 
ellipticity. 



We use an input model consisting of a lens of mass M = 7.5 x 
10"^ h-^MQ at redshift zi = 0.6 in a cosmological model with w = 
—0.8. We consider its critical curves for source redshifts z,s = 1 and 
Zs = 2, mimicking the constraints from two tangential arcs, and we 
fit their positions by varying the equation of state of dark energy 
and the lens mass. For simplicity, we consider only cosmologies 
with time-independent w. The fit is performed by minimising 



2 fxi{M,w)—xi 
Xi = 



X2(M,w) -X2 
'^2 



(22) 



where xi and X2 are the positions of the tangential critical curves 
of the input model for sources at Zs = 1 and Zs = 2, respectively, 
Al and A2 are their respective errors, and xi(M,w) and X2{M,w) 
are the corresponding positions of the critical curves predicted by 
the fitting model with mass M in a cosmological model with dark- 
energy equation of state w. We assume here to be in an idealized 
situation where the location of the critical curves is known at the 
1% level. 

We show the confidence levels in the w-M plane resulting from 
this fitting procedure in the upper left panel of Fig.|5| The inner- 
most and the outermost contours correspond to probability levels of 
68% and 99.7%, respectively. As anticipated, a good fit to the po- 
sition of the critical curves is obtained for a range of M and w, with 
99.7 confidence limits ranging between 6 x lO'"* /i^ ' Mq <M < 9 x 
lO^'^h-^Ma and -l<w< -0.65 (see also IChiba & Takahashil 



1200: 

For breaking this degeneracy, we must add constraints on the 
lens density profile. One possibility is to use the stellar veloc- 
ity dispersion data for the brightest cluster galaxies in conjunc- 
tion with the arc positions and redshifts (see Miralda-Escude 199^ 
ISand et al.''2003. and references therein). Including this additional 
constraint in our fitting procedure, we re-define our variable as 



3.4 Degeneracies in NFW models 

Since the expansion of the critical curves as a function of the source 
redshift is different for halos with different concentrations, a degen- 
eracy between mass and cosmology arises in NFW halos. That is, 
the same displacement of the critical curves for sources at two red- 
shifts can occur with halos of different concentrations in more than 
one cosmological models. As shown in Sect. 13. II the concentration 
is related to halo mass in NFW models of clusters. 

For probing this degeneracy, we carry out the following test. 



X2 — Xi + 



(23) 



where &;- is the velocity dispersion of the input model at lO/j^'kpc 
from the centre, A^^. is the uncertainty in its measurement, and 
ar{M,w) the value predicted by the fitting model. For calculating 
the velocity dispersion from the density profile we use the spherical 
Jeans equation, assuming isotropic orbits. 



dp(r)o.(r) GM{r)p{r) 



dr 



(24) 
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Figure 5. Confidence levels in the w-mass plane. The input model is a lens with an NFW profile with mass M = 7.5 x 10 Mr:, at z = 0.6 in a cosmological 
model with w = —0.8. The contours shown by the dark curve correspond to a probability level of 99.7%. The three panels correspond to constraints from two 
arcs at z = 1 and z = 2 (left); from two arcs at z = 1 and z = 2 plus velocity dispersion at 10 kpc from the centre (middle); from three arcs at z = 0.8, z = 2 and 
z = 3 plus velocity dispersion at \Qkpc from the centre (right). Note that these constraints apply only for the smooth analytical model we have used. 



where G is the gravitational constant and M(r) is the mass enclosed 
within a sphere of radius r. Again, we assume that the velocity 
dispersion at the chosen radius is known with an accuracy of 1%. 
This is a rather optimistic assumption, but our primary aim here 
is to discuss what kind of constraints are required for breaking the 
degeneracy. The new map of the confidence levels is shown in the 
upper right panel of Fig. |5| Although the contours shrink, some 
degeneracy remains, indicating that further constraints on the lens 
density profile are needed. 

Finally, we use a lensed image from a third source at red- 
shift Zs = 3 for locating another critical curve, in addition to two 
at Zs = 0.8 and Zs = 2 and to the velocity dispersion constraints. 
This allows us to distinguish among different cosmological models 
becomes easier as shown in the right panel of Fig.|5| but does re- 
quire that we get lucky with the arc redshifts. It is also valid only for 
the smooth mass distribution represented by our analytical model; 
real clusters are likely to have a more lumpy structure. 

Our simplified analysis with analytical NFW models for clus- 
ter halos was aimed at understanding the sensitivity of arc locations 
to variations in several physical parameters. Our results show that 
to get robust constraints on cosmology from a single cluster, a de- 
taile d modelling of the lensing mass distri bution is required (see 
also IChiba & Takaha shi 2002; Dalai et al.'2003'). This modelling 
would include the contribution from sub-structure in the cluster 
which we have not studie d in our analytical models. Further, as 
shown bv lDalal et alji200 4). structures along the line of sight pro- 
vide a source of error in modelling individual clusters that is very 
difficult to overcome. 



4 NUMERICAL MODELS 
4.1 Cluster sample 

Since asymmetries and substructures play a crucial role in 
determining the strong lensing properties of galaxy clusters 
iMeneghetti et aljl2003l) . analytic models can only be used for an 
approximate description of their lensing properties. More realistic 
mass distributions of clusters, as provided by numerical simula- 
tions, are needed for drawing quantitative conclusions. We now re- 



peat the analysis previously applied to analytic models to a sample 
of numerical clusters. 

The sample of clusters we use here consists of the 17 ha- 
los used by Dolae et al. (2004) and Meneahetti et al. (2004). We 
briefly summarise their main properties. The cluster models are 
pure dark-matter halos and were obtained using the most recen t ver- 
sion of the cosmological code GADGET Springel et alj200ll) . The 



code was extended bv lDolag etZI<2004r to cosmological models 
with dynamical dark energy. Each cluster was obtained by resimu- 
lating at higher resolut ion a patch of a pre- existing large-scale cos- 
mological simulation iTormen et alJll99'7l) . The initial conditions 
were set up with the purpose of obtaining identically comparable 
clusters in all the cosmologies at redshift z = 0. Of course, at higher 
redshifts, each of them appears at different evolutionary stages in 
different cosmologies, since the growth of the density perturbations 
depends on the equation of state of the dark energy. For example, 
clusters form earlier in the SUGRA and the DECDM compared to 
the RP and the ACDM models. The implications of the different 
time evolution of these clusters for different equations of state of 
dark energ y on the abundance of str ong lensing events is discussed 
in detail bv lMeneehetti et alji200 4h. 

The clusters in the sample contain on average Ny ~ 200, 000 
dark matter particles within their virial radii. The corresponding 
virial masses range between My = 3.1 x lO^'* to 1.7 x lO'^ h^^ Mq 
at redshift zero. 



4.2 Lensing simulations 

For each of them we take the snapshot at redshift z\ ~ 0.6 and per- 
form lensing simulations. We first select from the simulation box 
a cube of comoving side length 3/i^'Mpc. The three-dimensional 
mass distribution of the cluster is projected on a regular grid 
of 256 X 256 cells along three orthogonal axes. The resulting 
two-dimensional densit y fields are smoothed using t he Triangular 
Shaped Cloud method faockn ev & Eastwoodll98^ for avoiding 
discontinuities between neighbouring cells. Thus, we analyse 51 
projected lens mass distributions. 

Through each of the surface density maps, we trace a bundle 
of 1024 X 1024 light rays on a regular grid covering the central 
1 /i^'Mpc^ of the lens plane. This is large enough for enclosing the 
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Figure 6. Example showing how the critical curves of a numerical cluster at zi = 0.6 simulated in different cosmological models change between Zs = 1 and 



critical curves of all our numerical models. Deflection angles are 
computed using the method described in Meneahetti et al. (2000D. 
We first define a grid of 128 x 128 "test" rays, for each of which the 
deflection angle is calculated by directly summing the contributions 
from all cells on the surface density map S/.j, 



AG 



(25) 



where A is the area of one pixel on the surface density map and 
and Xij are the positions on the lens plane of the "test" ray (h, k) and 
of the surface density element Following IWambseanss et all 
il998h . we avoid the divergence when the distance between a light 
ray and the density grid-point is zero by shifting the "test" ray grid 
by half-cells in both directions with respect to the grid on which 
the surface density is given. We then determine the deflection an- 
gle of each of the 1024 x 1024 light rays by bi-cubic interpolation 
between the four nearest test rays. 
The reduced deflection angle is 



ax 



(26) 



Since we wish to estimate the growth of the tangential critical 
curves as a function of source redshift, we consider sources be- 
tween Zs = 1 and Zs = 2. 

Since the reduced deflection angle is the gradient of the lens- 
ing potential, 



V\(f(f ) = d{x) , 

the convergence and the shear can be easily written as: 



1 /dai da2 



(27) 



(28) 



Yi(^) = 2 
l2{x) = - 





da2\ 
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(30) 



The tangential critical curves are then determined by searching for 
those points in the lens plane where Eq. j9} is satisfied. 



4.3 Critical-curve sizes 

The critical curves for one of the clusters in our sample in the dif- 
ferent cosmological models are shown in Fig.|6| The sources are at 
redshift Zs = 1 in the upper panels and at Zs = 2 in the lower panels. 
Two important features are evident. First, for sources at redshift Zs 
the sizes of the critical curves differ substantially among the vari- 
ous cosmological models. For example, the cluster has almost no 
critical curves in the ACDM model, while they are already well de- 
veloped in the SUGRA model. The RP and the DECDM models fall 
between these two cosmologies. This is a consequence of the earlier 
formation epoch of clusters in the RP, DE CDM and in the SUGRA 
models than in the ACDM model (see e.g. | Bartelmann et alj2002L 



l2003tlDolag et all2004rMeneghetti et alJ2004 . due to which they 
have a larger concentration enabling them to be efficient lenses even 
at relatively high redshifts or for relatively close sources. Second, 
as expected from the analytical calculations, the relative enlarge- 
ment of the critical curves is higher in the ACDM than in the other 
cosmological models. 

We apply the same method used for the analytical models 
for estimating the size of critical curves of the numerical clusters. 
Fig. shows the relative growth of the critical curves in the four 
cosmological models as measured in the numerical simulations. 
Each curve represents the median among the 51 halos which de- 
velop a critical curve for source redshift Zs = 1. The number of 
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Figure 7. Critical curve sizes normalized to z, = 1 , as in Fig.|2] but for a nu- 
merical cluster sample comprising 17 halos placed at redshift z\ = 0.6. Only 
clusters producing critical curves for Zs = 1 are considered. Curves display 
the sample medians. The error bars show the first and the third quartile of 
the distribution, rescaled assuming that the information from 1000 pairs of 
arcs could be combined. 



useful clusters for this analysis ranges between ~ 20 in the ACDM 
to ~ 30 in the SUGRA model. The results confirm the qualitative 
expectations from analytic models: namely, the trend for different 
cosmologies. The absolute values of the relative growth are also 
consistent with the predictions for a moderate ellipticity (e<0.3) 
lensing potential (compare Figures 2 and 7). 

These results show that the statistical application of this 
method is potentially powerful. Upcoming s urveys from space, lik e 
those which will be conducted by SNAP ( Alderine et a i] |2004) . 
could provide detailed observations of order thousand galaxy clus- 
ters, allowing the information from many lenses to be combined. 
The error bars in Fig. Q show the first and the third quartiles of 
the curve distribution we obtain from our numerical cluster sam- 
ple. They were rescaled to the expected error when the information 
from ~ 1000 pairs of arcs is combined. The figure shows that when 
constraints on the position of the critical curves from sources at 
significantly different redshifts and in a sufficiently large sample of 
clusters are used, it becomes possible to discriminate among differ- 
ent cosmological models. We have used a simple choice of source 
redshifts; by extending the analysis to a wider redshift range, es- 
pecially to redshifts beyond 2 for the distant arcs, the constraints 
became stronger. The analysis needs to be extended in other ways 
as well, by combining the information from different lens redshifts 
and finding the best way to weight a given cluster. 

In our test, we have used a simple estimate of the size of 
the critical curves, which we can assume is related to the cluster- 
centric distance of an observed arc, without any mass modelling 
of the cluster. We obtain substantially different amounts of growth 
of the critical curves depending on the level of substructure in the 
central region of the clusters. If the cluster is close to critical for 
sources at redshift ~ 1 at the location of a substructure, we mea- 
sure a rapid growth of the critical curve once the sources are shifted 



to higher redshift. This occurs because the critical surface density 
becomes smaller and the cluster critical curves wiggle around the 
emerging critial mass lumps. This is shown for example in Fig.|5| 
where the shape of the critical curves changes dramatically between 
Zs = 1 and Zs = 2 for some of the cosmological models (see e.g. the 
DECDM model). If mass modelling for some of the clusters were to 
be included, the scatter in Fig.Qcould be reduced. Finally we note 
that to some extent carrying out such an analysis from observations 
would require simulations that correctly reproduce the properties of 
real clusters. While we have used only the relative locations of arcs 
at different redshifts, it is worth examining what aspects of cluster 
structure these are sensitive to. 



4.4 Imaging surveys for cluster lensing 

Strong lensing studies with galaxy clusters require high quality 
imaging data. While we will not examine the survey requirements 
in any detail, we have used our numerical clusters to produce lensed 
images of the Hubble Ultra Deep Field (HUDF). We can use these 
to study of the effects of sky brightness, photon noise and resolu- 
tion on cluster arcs. 

Figure 8 shows the kind of imaging expected from space imag- 
ing and current state of the art ground based imaging. Clearly, the 
demands on resolution are very high to reliably measure the prop- 
erties of lensed arcs. With 0.5 arcsecond imaging the quality of 
the images is compromised, as shown in the right panel, but the 
large arcs are still identified. Adaptive optics may enable better res- 
olution images to be obtained than shown in this figure; short of 
that, it may still be possible to use giant arcs and on average get 
the arc separations accurately enough. Ideally, a space based sur- 
vey that images a large number of multiple-arc clusters with reso- 
lution comparable to the HST would provide an adequate sample. 
For such a survey, the question of how deep it is worth going to find 
high-redshift arcs merits worth further study. 

In addition to multi-color imaging, which would enable photo- 
metric redshift estimates, cluster cosmography would require spec- 
troscopic redshifts for the lensed arcs. An important exercise for fu- 
ture work is to estimate the number of arcs expe cted from a realistic 
sample of clusters for given survey parameters. iDalal et alj J2003h 
have used three existing cluster lensing datasets ( Gladde rs et alj 
2003; LuDDinoetal. 1999; Zaritskv & Gonzalez 2003) and esti- 
mated that of order 1000 arcs are expected on the sky even from 
a survey of modest depth. Given the landscape of proposed tele- 
scopes and surveys, it is not prohibitive to think high resolution 
imaging of order a thousand cluster containing fields, as well as 
follow-up spectroscopy on the several thousand arcs that might be 
found. 



5 SUMMARY 

This paper presents a study of cosmography with galaxy clusters 
that produce lensed arcs of background galaxies. Clusters with mul- 
tiple arcs from different redshifts provide constraints on angular di- 
ameter distances at different redshifts. These in turn constrain the 
equation of state of the dark energy. We have studied this effect 
with analytical models to understand some of the qualitative trends, 
and with numerical simulations to see what can be expected from 
surveys of real clusters. 

We have used analytical models of clusters based on the NFW 
profile to study how the positions of arcs depend on the mass 




Figure 8. Simulated observations from a VLT-like telescope of the same field observed by HST in the ultra-deep field with 10* sec exposure. The simulations 
include a sky brightness / = 19.9 mag/arcsec", Poisson photon noise and a seeing of 0.5 arcsec. The lensing cluster of galaxies is at = 0.25. The cluster 
members have been added by placing a cD galaxy at the cluster center and by distributing the other galaxies following the underlying dark matter mass 
distribution. The galaxy luminosity reflects that of realistic cluster members (spectral energy distribution, morphological type as a function of the the cluster- 
centric distance reproduce observational constraints). The integration time is 3000 sec. In the lensing simulation, all the sources were assumed to be at ^ = 2. 
The field size is 1 arcminute along the x-axis. The middle panel includes sky brightness and Poisson photon noise; the right panel includes seeing. 



distribution and the cosmological model. We find that the ellip- 
ticity and concentration of the cluster halo can shift the criti- 
cal curves. Substructure would have the same effect. This means 
that for individual clusters, it is important to have independent 
information on the mass distribution, especially inside the criti- 
cal curve, to extract cosmographic information from arc positions 
and redshifts. Our analysis overlaps with o t her recent studies o f 
cluster arcs, e.g. £hiba & Takahashj f2002h: lOguri et alj l200?t): 
iMeneghe tti et al j 12003.) : ■Dalai et all ilOOshn^^mbsgans^^l] 
(2004). We explore examples of how information from inner ve- 
locity dispersions and multiple arcs could be combined for a single 
cluster. However a more sophisticated study is warranted. 

Given the limitations of analytical models (they cannot for ex- 
ample include the effect of substructure in the halos), we conclude 
that numerical modelling is essential to extract information from 
observed clusters. We use a sample of simulated clusters in four 
dark energy models in an exploratory study of whether future sur- 
veys can provide cosmographic information. We examine the scat- 
ter in the relative positions of arcs due to differences in the mass 
distributions of individual clusters. We focus on the critical curves 
and use these to estimate arc positions. While this is a simplified 
study, we believe it gives us a reasonable estimate of the minimum 
requirements of a cluster sample to overcome the cluster-to-cluster 
scatter and obtain cosmographic constraints. 

In the numerical study we have not attempted to model the 
mass distribution of individual clusters. We have only used infor- 
mation on the relative sizes of critical curves at different redshifts. 
In practice, some constraints on the mass distributions will be ob- 
tained from strong and weak lensing, as well as X-Ray, SZ and 
velocity dispersion measurement. Further study is needed on how 
well this information can be used to improve cosmographic con- 
straints. Further, we have used information only from arcs up to 
redshift 2; one would do better if imaging and spectroscopy on arcs 
at higher redshifts were available. 

We have attempted to include the main sources of complexity 
that arise from the mass distribution of clusters. We have neglected 
several factors in the observation of clusters and their interpreta- 
tion that merit further study. Structures along the line of sight at 



redshifts different from the cluster are an important source of scat- 
ter. Dalai et al. 1 2004) have studied the impact of these structures; 
they find that they limit the use of individual clusters and could 
contribute to the error in dark energy parameters from an ensemble 
of clusters. Further we used just the positions of the largest arcs, 
and have largely ignored other strong lensing information as well 
as the measurement error in using arc positions. These factors are 
all important and require more detailed studies. 

We have estimated the scatter in cosmographic constraints 
from different clusters. We have not addressed the question of how 
an actual survey would best average results from different clusters. 
Clearly some clusters will have multiple arcs (over 100 in the case 
of Abell 1689) and therefore contain more information. Others may 
have a regular, spherically symmetric structure which would make 
it easier to extract cosmographic information from arc positions. 

While completing this study, we have learnt of a preprint 
by Dalai et al. (.20041) that examines sources of noise in arc- 
cosmography. 
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